The Strange Behavior of Critical Branched Polymers 
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We find that 2-dimensional (2-D) critical branched poly- 
mers with no impurities conclusively belong to the same uni- 
versality class as 2-D random percolation clusters, although 
pure critical 3-D branched polymers do not belong to the 
3-D percolation universality class. We find, moreover, that 
the fractal dimension df of critical branched polymers is in- 
dependent of the presence of a random environment in 2- 
D, but not in 3-D. We also report that when there are 
no impurities the critical branching probability in 3-D is 
6c = 3.34 X 10~* ± 0.16 X 10^*. 

Polymerization has been an important topic of recent 
interest in the physics community Polymers are 

extensively studied ||^-^ complex systems, and can be 
either linear or branched p],p|-p^ . In spite of being sim- 
pler, linear polymers (LP) present physical features such 
as non-Markovian growth and nontrivial scaling. For sev- 
eral decades it has been known that linear polymers in 
diluted solutions can be successfully modeled using self- 
avoiding walks (SAW). More recently it was found that 
large SAW chains can be more easily constructed by ki- 
netic growth walks (KGW) [0 that belong to the same 
universality class as SAW chains [Q. Branched poly- 
mers (BP) have a larger degree of complexity and exhibit 
an astonishingly rich phenomenology ||l|,^,^,|6P, [lO| , [l^ -p^ . 
In the 1990's, the KGW approach was generalized to de- 
scribe the growth of BP in disordered media |p], p^ , p^ . 
Such BP are typically grown when polymerization occurs 
in a solution with two types of units: (i) monomers with 
two chemically active molecular sites ("tips") that lead 
to linear polymerization and (ii) three-tipped mononers 
that lead to branchings (i.e., bifurcations). Many impor- 
tant physical properties of the resulting BP depend on 
the probability b for branching to occur and the concen- 
tration c of impurities that can block or slow the 
polymerization. 

Here, we first address the controversial question: to 
which universality class does the critical polymer (6 — be) 
belong for c = 0? A typical phase diagram p|, p^ , p^ 
(Fig. 1(a)) for BP consists of several distinct regions. For 
large b and small c, the polymer is compact and grows 
indefinitely (i.e., is "infinite") with a smooth (faceted) 
growth surface, while for smaller values 1 > 6 > 6c closer 
to the critical line 5c (c) the resulting polymer grows in- 
definitely but with a rough growth surface ||l6|. For yet 
smaller values of 6 < 5c the resulting polymer is finite. 



due to termination of polymerization at all active growth 
sites. Along the critical line b — bc{c), the polymer is 
fractal (Fig. 1(a)) and has a diverging correlation length. 

The commonly held belief that BP belong to the lat- 
tice animals [Q universality class was based on the as- 
sumption that there is no random environment (c — 
0) and that there are only repulsive forces between 
chains 1151] . It was thus a remarkable finding when Bunde 
et al. |l| showed, using the branched polymer growth 
model (BPGM), that critical BP in random environments 
in 2-D belong to the same universality class as 2-D ran- 
dom percolation clusters. One conceivable explanation 
for this interesting result is that in 2-D such behavior is 
induced by the random impurity obstacles that, indeed, 
form real percolation clusters with site occupation prob- 
ability p ~ c. Hence, it can be argued that the polymer 
grows as if constrained by the random obstacles, and this 
effect can generate an effective attraction between chains 
(Fig. 1(b)). For the special case c — 0, however, there 
are no such impurities, hence this physical argument be- 
comes irrelevant. The study by Bunde et al. is conclu- 
sive all along the critical line b = 5c (c) for c 7^ 0, but the 
important special case of pure BP (c = 0) still remains 
unconfirmed |l|, and both the percolation as well as the 
lattice animals universality classes remain plausible. An 
additional controversy surrounding the special case c — 
arises because in many systems p5| the introduction of 
disorder is sufficient to modify the universality class (and 
sometimes even to destroy the ordered phase) . The Har- 
ris criterion predicts that the critical BP in 2-D is a 
candidate for experiencing a change of universality class 
when going from c = to c ^ 0. 

We investigate this problem by using simulations of 
the BPGM 0. BPGM is a generalization of KGW |l|] 
which is able to capture the essential dynamics of 
branched polymerization. We briefly describe the model 
here (Fig. 1(b)). The BP is generated from an initiating 
seed at time t = located at the center of a d-dimensional 
lattice. At time step t = 1, one of the vacant nearest- 
neighbor sites of the seed is chosen randomly and occu- 
pied to become the next active growth tip, and the seed 
ceases to be an active growth tip. At each step t of the 
growth process, a branching also can occur at every such 
growth tip, by occupation of an additional (second) site 
with probability 5. A fraction c of randomly chosen sites 
is not available for growth. Hence, at time t-i-1, the poly- 
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mer can grow from any of the active tips added at the 
previous step t to empty neighboring sites (Fig. 1(b)). 

Our extensive 2-D simulations, discussed below, leave 
no doubt that, surprisingly, random percolation is indeed 
the correct universality class even for c — 0. We simu- 
late critical polymers for c = and b = bc{0) « 0.056 
and find that the fractal dimension d{, the minimum di- 
mension dmin (Fig- 2(a)) and the chemical dimension di 
(Fig. 2(b)) are in strong agreement with the known values 
for a a 2-D critical percolation cluster flj] . These dimen- 
sions are defined by the relations M ^ r^f , M ~ i"^' , and 
£ ^ T-'imin ^ -^^here r is radius and £ is the chemical distance 
(£ = tm BPGM). Clearly di = ■ d^in ~ 1-89. We con- 
clude that the critical 2-D random percolation cluster and 
the critical polymer generated by BPGM unmistakably 
belong to the same universality class. Therefore the frac- 
tal dimension of the critical BP in 2-D is independent 
of the presence of a random environment — in apparent 
violation of the Harris criterion. 

We further investigate this "paradoxical" result by per- 
forming similar 3-D simulations of pure (c = 0) critical 
BP, because the Harris criterion favors a little more the 
irrelevance of the quenched disorder in 3-D. Until now, 
such 3-D simulations have not been possible because the 
precise value of &c(0) in 3-D was unknown. 

So we now determine the value of &c(0) in 3-D for 
c — 0. The fundamental difficulty of finding a precise 
value for 6c (0) in 3-D is that the value is extremely small, 
&c(0) < 10~^ (i.e., on average more than 10^ monomer 
units between successive branchings), making it compu- 
tationally expensive to pinpoint it using traditional nu- 
merical approaches that involve tuning b. We overcome 
this difficulty by mapping BPGM to the computationally 
less intensive fixed number of tips model (FNTM) 
that has been shown llTUlSfl to generate phase diagrams 
identical to BPGM. FNTM differs from BPGM as fol- 
lows: rather than fixing the branching probability 6, in- 
stead FNTM dynamically attempts to fix the number N 
of active polymerization growth tips. Hence, the branch- 
ing process does not occur with an a priori fixed value 
b, but rather branchings occur whenever one or more of 
the N existing active tips "die," either because of impu- 
rities or due to other steric hindrance effects. We can be 
certain that the value of 6c found using FNTM is iden- 
tical to the one found from BPGM simulations because 
FNTM generates polymers in which the number of active 
growth tips neither vanishes nor explodes exponentially, 
i.e., corresponding to the critical line at 6 = 6c in BPGM 
simulations |17 lq|. FNTM greatly reduces the computa- 



tional burden of generating critical polymers because it 
automatically generates critical polymers with 6 = 6c (c) 
and, moreover, it has only one free parameter, namely 
the impurity concentration c. (The arbitrary value of 
iV > 1 chosen is irrelevant in the large mass limit, as 
seen in Fig. 3(a)) The effective critical branching proba- 
bility 6cff = B/M is found by dividing the total number 



B of branchings of active tips by the total number M of 
occupied polymer sites. Most importantly, FNTM goes 
spontaneously to the critical line without the need of pa- 
rameter tuning (of 6), as with BP growth models that 
exhibit self-organized criticality (SOC) ]l8[ |. 

We simulate FNTM for c = in 3-D and N = 
5, 10, 20, 30, 50, 80 and 100 tips as a function of 
polymer mass M. In the limit M — > oo, we find that 
the values of 6cff = B/M converge towards an identi- 
cal point independently of N (Fig. 3(a)). The uncer- 
tainty for a typical point in Fig. 3(a) is A6eff « 10~^. 
This remarkably small error bar is made possible only 
because FNTM dynamically seeks the fixed point attrac- 
tor near 6 « 6c, such that the effective critical branch- 
ing probability 6cff — B/M fluctuates around 6c. By 
extrapolating b^s for 1/M ^ we flnd that 6c(0) = 
3.34 X 10-'i±0.16 X 10-''. 

Since it is not altogether inconceivable that the value 
of 6c (0) obtained from FNTM will not carry over to 
BPGM, we check the above conclusions by simulating 
BPGM with 2.9 X IQ-* < 6 < 4.0 x 10"''. We find 
(Fig. 3(b)) that the number of active tips either grows 
or decays much faster than a power law except in the 
approximate range 3.3 x 10""* < 6c(0) < 3.4 x 10""*. This 
value obtained for 6c (0) is not inconsistent with the value 
6c = 3.34 X 10-'* ± 0.16 X 10-°-* found using FNTM. 

Using this result, we are able to perform 3-D BPGM 
simulations of pure critical BP. We find that d{ w 
dmin ~ 2 and w 1 in 3-D for 6 = 6c (0) and c = 
(Figs. 2(c), (d)), definitely ruling out the 3-D random 
percolation universality class. Although this is the ex- 
pected field theoretical result (see ref. ||l^) for 3-D poly- 
mers, the Harris criterion indicates a relevant disorder 
and a change in the value of dj (because diy = 1.5 < 2 
for d — 3 BP). This prediction is confirmed by the finding 
of Bunde et al. jl^, that d{ « 2.53 for 6c — 1, consistent 
with 3-D random percolation. 

Somewhat surprising is our confirmation that in 2-D 
the fractal dimensions of the critical BP are identical for 
the pure case c = and the maximally disordered case 
6 = 1, c « 1 — Pc (see Fig. 1(a)), where Pc « 0.59 is the 
critical site occupation probability in random percolation 
for a square lattice. For 6=1 the critical polymer grows 
in a vacant space that is itself almost a critical percola- 
tion cluster (Fig. 1(a)). Since there is no random envi- 
ronment for c = 0, a larger fractal dimension had been 
conceivable in principle. Our findings indicate, strangely, 
that the fractal dimension (but not the lacunarity) of the 
critical polymer in 2-D is independent of the growth envi- 
ronment. Why this is so is indeed an interesting question 
that merits further investigation, and seems to be related 
to a self-organizing process Criticality occurs only 

when there is a delicate balance between the rates of tip 
deaths and branchings — a dynamic not unlike the one 
used in the self-organized generation of infinite 2-D per- 
colation clusters 19 1 . 
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We now comment briefly on our finding that, strangely, 
pure critical BP fall into the percolation universality class 
in 2-D but not in 3-D. The deeper reason behind this mys- 
terious behavior is possibly related to the well known fact 
that in 2-D, encounters between different branches of a 
BP can result either in "trapping" of a polymer chain 
or else in "scattering" (bending) of the chains. In 3-D 
a completely different effect becomes important in the 
large-scale growth dynamics: interpenetration and "en- 
tanglement" of BP chains j|] — not by independent linear 
chains, but by linear parts of a BP that act independently 
at criticality, because the chemical distance between suc- 
cessive branchings is much larger than the persistence 
length. In 2-D, there is no entanglement (for topologi- 
cal reasons), but in 3-D such effects can counteract the 
excluded volume repulsion, so that the mean square dis- 
placement of the tips grows linearly with mass (as with 
polymer melts, see Ch. 2 and 8 of ref. (D), leading to 
power law scaling equivalent to ideal Gaussian chains 
with df — 2. Effectively, there is a larger repulsion be- 
tween chains in 2-D, thereby leading to superdiffusive 
motion of the tips (with d{ < 2). Note that trapping in 
3-D is extremely rare, hence the very low value 6c (0). 

Finally, one hint for the cause of the observed change 
in the universality class (when we vary b from b = 6c (0) 
to 6 = 1) in 3-D becomes evident by comparing the two 
situations with zero and maximum disorder: for 6c (0) = 
3.3 X 10~^ the "freedom" of motion of the tips is limited 
only by the sparse steric hindrance due to the polymer 
chains, while for 6 = 1 the behavior is dictated by the 
intersticial percolation cluster available for BP growth. 
The arguments presented above raise new questions, such 
as whether other critical exponents vary along the critical 
line 6 = 6c (c). 

We thank U. L. Fulco, H. J. Hilhorst, M. L. Lyra, 
S. Roux, and L. R. da Silva for discussions; CNPq, 
PRONEX, the Projeto Nordeste de Pesquisa & FINEP 
for support. 



[1] A. Bunde, S. Havlin and M. Porto, Phys. Rev. Lett. 74 
2714 (1995). 

[2] P. De los Rios, Phys. Rev. Lett. 82 4236 (1999). 

[3] M. Myllys at al., Phys. Rev. Lett. 84 1946 (2000). 

[4] R. Du et al., Phys. Rev. Lett. 84 2417 (2000). 

[5] P. J. Flory, Principles of Polymer Chemistry (Cornell 
University Press, Ithaca, 1953). 

[6] P. G. de Gennes, Scaling Concepts in Polymer Phystcs 
(Cornell University Press, Ithaca, 1979). 

[7] M. Doi and S. F. Edwards, The Theory of Polymer Dy- 
namics (Oxford University Press, Oxford, 1986). 

[8] T. C. Lubensky and J. Isaacson, Phys. Rev. Lett. 41, 829 
(1978). 

[9] G. Parisi and N. Sourlas, Phys. Rev. Lett. 46, 871 (1981). 
[10] L. S. Lucena, J. M. Araiijo, D. M. Tavares, L. R. da Silva 
and C. Tsallis, Phys. Rev. Lett. 72, 230 (1994). 



[11] 

[12] 
[13] 

[14] 

[15] 

[16] 
[17] 



[18] 
[19] 



I. Majid, N. Jan, A. Coniglio and H. E. Stanley, Phys. 
Rev. Lett. 52 1257 (1984). 

L. Peliti, J. Phys. (Paris) Lett. 45 L925 (1984). 

M. Porto, A. Shehter, A. Bunde and S. Havlin, Phys. 

Rev. E 54 1742 (1996). 

A. Bunde and S. Havlin, Fractals and Disordered Systems 
(Springer, Berlin, 1991). 

J. Cardy, Scaling and Renormalization in Statistical 
Physics (Cambridge University Press, Cambridge, 1997), 
and references therein. 

L. S. Lucena, L. R. da Silva and S. Roux, Physica A 266 
81 (1999). 

H. H. Aragao-Rego et al., in Computer Simulation Stud- 
ies in Condensed Matter Physics XII , eds. D. P. Landau, 
S. P. Lewis and H.-B. Schiittler, 85 223 (Springer, Berlin, 
1999). 

J. S. Andrade Jr., L. S. Lucena, A. M. Alencar and J. E. 
Freitas, Physica A 238 163 (1997). 

A. M. Alencar, J. S. Andrade Jr. and L. S. Lucena, Phys. 
Rev. E 56, R2379 (1997). 



bjc) fractal 

v. ^ 




(b) 



FIG. 1. (a) schematic phase diagram for BP (see text). 
We focus exclusively on the critical point at 6 = fec(O) that 
corresponds to BP growth in the realistic case of no impu- 
rities, (b) the branched polymer growth model (BPGM) as- 
sumes that growth takes place only at active growth tip sites 
located on a lattice. The seed (cross) initiates polymeriza- 
tion governed by a branching probability b. When the im- 
purity concentration c is nonzero, then the random obstacles 
(circles) can sometimes cause an effective attraction between 
chains (see branches on the right), but this effect disappears 
when c = 0. 
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FIG. 2. (a) double-log plot of the polymer mass in units 
of monomers (mass axis label at left) and chemical distance 
(label at right) as a function of radius for BP with critical 
branching probability 6c (0) in the absence of impurities in 
2-D. We have used 100 samples with mass Af = 5 x 10^ 
monomers. The fractal dimension df and the minimum di- 
mension dmin (the slopes) are those of a 2-D critical per- 
colation cluster (straight lines), (b) double- log plot of the 
polymer mass as a function of the chemical distance for the 
identical simulations, showing that the chemical dimension is 
also consistent with the universality class of 2-D random per- 
colation, (c) polymer mass and chemical distance vs. radius 
and (d) polymer mass vs. chemical distance in 3-D, show- 
ing that the behavior is definitely different from 3-D random 
percolation (df = 2.52, di = 1.84, dmin = 1-37 @). We 
have averaged over 900 samples. We find that 3-D critical BP 
scaling is more reminiscent of polymer melts Q than random 
percolation. 
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FIG. 3. (a) fecff a function of 1/M and the number A'' 
of active growth tips for FNTM in 3-D. The solid lines are 
fitted curves, by which we are able to extrapolate for arbi- 
trarily large masses. The values shown for 1/M = are the 
extrapolated results. We are thus able to estimate the value 
of 6c (0) = 3.34 X 10~* ±0.16 x 10"". (b) double log plot 
of A?^ as a function of chemical distance I for BPGM, with 
A^o = 2 tips initially at t = 0. For values 6 larger than the 
estimated value of 6c (0), we find that A'' tends to explode 
exponentially, while for 6 < 6c (0) we find TV decays. For 
3.3 X 10"" < 6 < 3.4 X 10"", however, we obtain power law 
growth for A', indicating critical branching. 
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